FINITE HEAT KERNEL EXPANSIONS ON THE REAL 

LINE 
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s ! ' Abstract. Let £ = d 2 /dx 2 + u{x) be the one-dimensional Schrodinger 

operator and H(x, y, t) be the corresponding heat kernel. We prove that 
■ the nth Hadamard's coefficient H n (x,y) is equal to if and only if there 

exists a differential operator M of order 2n — 1 such that £ 2n_1 = Jvl 2 . 
Thus, the heat expansion is finite if and only if the potential u(x) is a 
00 ' rational solution of the KdV hierarchy decaying at infinity studied in 

DHJ. Equivalents, one can characterize the corresponding operators C 
as the rank one bispectral family in [Sj. 

> 

O ■ 1- Introduction and description of the main result 

0> Let u{x) be a smooth function and let C be the one-dimensional Schrodinger 

operator 

<? 



in 
o 



c = d^ +u{x) - (L1) 

The heat kernel H(x,y,t) is the fundamental solution of the heat equation 



22 

(l-C)f = 0. (1.2) 



X 



at 

It is well known that H(x, y, t) has an asymptotic expansion of the form 



H( X , y, t) ~ + H ^ x ' y^ n J as f -* 0+ ' ( L3 ) 

valid for x close to y. The coefficients H n are the so called Hadamard's 
coefficients [IT] . 

In a series of papers Berest and Veselov studied linear hyperbolic partial 
differential operators that satisfy the Huygens principle and discovered a 
beautiful connection with integrable systems and the bispectral problem jS], 
see for example jlj |3] and the references therein. Grunbaum jH| observed 
that a similar connection is present in the context of the heat equation. More 
precisely, he showed that for the first few potentials = 2 log(rj(x))", 

with ro = 1, T\{x) = x, T2(x) = x 3 /3 — S3, formed by the rational solutions of 
Korteweg-de Vries (KdV) equation, the asymptotic expansion gives rise 
to an exact formula consisting of finite number of terms. The corresponding 
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operators C belong to the rank one family of solutions to the bispectral 
problem 

The main result of the present paper is to prove that this holds for all rank 
one bispectral operators and that this completely characterizes the rank one 
family. More precisely we have 

Theorem 1.1. The following conditions are equivalent 

(i) H n (x, y) = for all x and y; 

(ii) There exists a differential operator M. of order 2n — 1 such that 

C 2n-1 =M 2 (L4) 

Operators C and M. , satisfying (|1.4|) were studied in a paper by Burchnall 
and Chaundy 1 6 . Their main result essentially means that the operator C 
can be obtained by applying a sequence of rational Darboux transformations 
to the differential operator Cq = d 2 /dx 2 . According to a result by Adler 
and Moser the corresponding potentials u(x) are precisely the rational 
solutions of the KdV equation decaying at oo. The latter were discovered 
by Airault, McKean and Moser |5j and mysteriously connected with the 
Calogero-Moser system. Finally, the operators C in (ii) coincide with the 
rank one bispectral family in the work of Duistermaat and Griinbaum (HJ. 
Thus, the present paper adds one more characterization of the operators C 
satisfying condition (ii), namely the finiteness property of the heat kernel 
expansion. 

The proof of this theorem can be obtained using the connection between 
the heat kernel and the KdV hierarchy. This goes back to the pioneer work 
of McKean and van Moerbeke ^Zj, where it was shown that, restricted on 
the diagonal x = y, the Hadamard's coefficients give the flows of the KdV 
hierarchy. See also [2j and the references therein for extensions to matrix- 
valued heat kernel expansions and applications. 

In a recent paper Jl] we showed that one can go also in the opposite 
direction and construct the heat kernel using the r-function (in the sense of 
Sato see |18(I13|). We combine this formula for H n (x, y) in terms of the Baker 
function, the explicit description of the coordinate ring for the algebraic 
curves and the bilinear (Hirota) equations to show that the vanishing of 
the Hadamard's coefficient H n (x,y) is equivalent to condition (ii) in the 
theorem. 

For related results concerning the finiteness property of the heat kernel 
expansion on the integers and rational solutions of the Toda lattice hierarchy 
see jlUj . For solitons of the Toda lattice and purely discrete versions of the 
heat kernel see |12| . 



They considered even a more general situation of operators P and Q of arbitrary 
orders, satisfying P dc s(Q) = Qd cg (P) 
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2. Baker functions and an infinite Grassmannian 

In this section we collect some preliminary facts about the Segal- Wilson 
Grassmannian |2f)j and the parametrization of rank one commutative rings 
of differential operators. We also present the formula from 14^ which relates 
H n (x, y) with Sato theory for the KdV hierarchy. Finally we prove a simple 
lemma which allows us later to use the algebro-geometrical data. 

Let us denote by 5 1 the unit circle S' 1 = {2:GC:|z| = l} and let H 
denote the Hilbert space L 2 (5 1 , C). We split H as the orthogonal direct sum 
H = H + © where H + and (resp. H-) consists of the functions whose 
Fourier series involves only nonnegative (resp. only negative) powers of z. 
We denote by Gr the Segal- Wilson Grassmannian of all closed subspaces W 
of H such that 

• the projection W — ► H + is a Fredholm operator of index zero; 

• the projection W — > H_ is a compact operator. 

Let Gr^ 2 ) be the subspace of Gr given by 

Gr (2) = {W G Gr : z 2 W C W}. 

To each space W G Gr there is a unique (Baker) function ^w(x, z) charac- 
terized by the following properties: 

• V w has the form V w {x, z) = (l + E£i M^z^) e xz ; 

• ^w(x, .) belongs to W for each x. 

We denote by W al& the subspace of elements of finite order, i.e. the 
elements of the form Y^j<k a 3 z3 • This is a dense subspace of W. A very im- 
portant object connected to the plane W is the ring Ayy of analytic functions 
f(z) on S 1 that leave W alg invariant, i.e. 

A w = {/(z) analytic on S 1 : f(z)W als C W alg }. 

For each f{z) £ A\y one can show that there exists a unique differential 
operator Lf{x,d/dx) such that 

L f W(x,z) = f(z)^ w (x,z). 

The order of the operator Lf is equal to the order of /. In general, Ayy is 
trivial, i.e. Aw = C. The spaces W which arise from algebro-geometrical 
data (which are important for us) are precisely those such that A\y con- 
tains an element of each sufficiently large order. The commutative ring of 
differential operators 

% w = {L f :f€A w } 

is isomorphic to Ay/. Every rank one commutative ring 2 of differential 
operators can be obtained by this construction for an appropriate W 6 Gr. 
Notice that if W € Gr (2) then z 2 G W and C = L z 2 = d?/dx 2 + u(x) is a 
second-order operator. 



A commutative ring of differential operators is called rank one if it contains an operator 
of every sufficiently large order. 
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Remark 2.1. In Sato's theory of the Kadomtsev-Petviashvili hierarchy one 
usually writes the Baker (wave) function 

depending on all times s\ = x, S2, S3, . . . . Then if W £ Gi^ the operator 
C = L z 2 satisfies the KdV hierarchy 

— = [(£ k /*) + ,£], k = 1,3,5,..., (2.1) 

where (£ fc//2 )+ is the differential part of the pseudo-differential operator C k / 2 . 

On the Hilbert space H we have the non-degenerate continuous skew- 
symmetric bilinear form 

{f,g) = ^-j J(z)g{-z)dz. (2.2) 

If W E Gr then its annihilator 

W* = {f e H : (f,g) = for all g € W} 
also belongs to Gr. The adjoint Baker function fy*(x,z) is defined by 

It follows immediately from the definition that the following bilinear identity 

/ * w (x,z)** w (y,z)dz = (2.3) 
Js 1 

holds (see [7]). 

The main result of |14j is the following formula for H n (x,y) 

H n (x,y) = 1 , 2w _ 1 / g n {2{x-y)z)^ w {x,z)^* w (y,z)dz, (2.4) 

iri(x - y) 2n 1 J s i 

where 



g n (z) = (-1) 



n-1 



t4 (2ra-fc-2)! fe 

^ kUn-k- l)\ Z 
.k=0 y ' 



e"5. (2.5) 



It is interesting to notice that the function g n (z) is closely related to a 
specific Laguerre polynomial. Recall that the Laguerre polynomials L^(z) 
are defined by the following formula 

L«(z)-^ T{n + a + 1) { ~ z)k (2 6) 

L ^ z )- 2^ r{k + a + 1) k \(n-k)V ^ 

where T(z) is the usual gamma function, see for example l i6, p. 77]. Thus 
we can write g n (z) as 

g n {z) = (n-l)\L- 2n + 1 (z)e-i. (2.7) 
Using the above formula we prove the following lemma. 
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Lemma 2.2. Let 

oo 

g n (z) = J2(3 k z k , (2.8) 

fc=0 

be the Taylor expansion of the function g n (z) around z = 0. Then fyj-i = 
for 1 < j < n 

Proof. It is well known that the function f(z) = L"(z)e~2 satisfies the 
differential equation 



2 r 4 



see for example ^3 p. 85]. 

Using the above equation and H2.7|) we see that g n (z) must satisfy the 
differential equation 

zg n -2(n-l)g n = — , 
from which it follows immediately that 

8(2j - - n)(3 2j -i = foj-s- (2.9) 
Combining this with ft = we obtain the assertion of the lemma. □ 

3. Proof of the theorem 

The implication (i)=^(ii) can be proved using the connection between 
Hadamard's coefficients on the diagonal and the KdV hierarchy ()2.1|) . In- 
deed we have 

8H k (x,x) _ 2 k ~ l r (c ^L. f] . , 

dx ~ (2k - l)!! i(£ )+ ' £J ' (3 - 1} 
see for example ^2 Theorem 5.2]. If H n (x, y) = 0, then we have Hk(x, y) = 
for all k > n. Formula 1)3.1 j) combined with (|2.1jl shows that u is stationary 
under the flows d/ds2k-i for k > n. This implies that Aw D C[z 2 , z 2 '" 1 ^ 1 ]. 
Thus the operator J\A = L z i n -\ satisfies M 2 = C 2n ~ l , which gives (ii). 

Below we prove the implication (ii) =>• (i). If (ii) holds then the operator C 
commutes with the operator M of order 2n — 1, i.e. it belongs to a rank one 
commutative ring of differential operators f° r some plane W £ Gi {2) . 
Moreover we have 

C(x, d/dx)^ w (x, z) = z 2 ^ w {x, z) (3.2) 
M(x, d/dx)^ w (x, z) = z 2n ~H w {x, z), (3.3) 

or equivalently 

C^z^C^. (3.4) 
This combined with 1)2. 3|) means that 

k is even or 
k > 2n - 1 



z^ w (x, z)^* w (y, z)dz = if 

s 1 



On the other hand, Lemma f2.2l savs that the first odd power of z in the Taylor 
expansion of g n (z) is z 2n ~ l . Thus formula 1)2.4)1 shows that H n (x,y) = 0. 
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Remark 3.1. The subspaces W parametrizing the operators C described in 
Theorem ll.lf ii) form a sub-Grassmannian of Gr^ 2 ) denoted by Gr^ in |2Dj . 
It also parametrizes all rank one commutative rings of differential operators 
21 with unicursal spectral curve SpecSl, containing an operator of order two. 
It is well known (see [201 p. 46]) that it has a cell decomposition with cells 
indexed by the sets 

S n = {—n, —n + 2, — n + 4, . . . ,n, n + 1, n + 2, . . . }. 

(2) 

The corresponding cell C n in Gi a consists of the solutions to the KdV 
hierarchy (|2.1|) flowing out of the initial value 

n(n + 1) 



u(x, 0, 0, ... ) 



,2 



These solutions exhaust all rational solutions to the KdV hierarchy that 
vanish at x = oo ^H]- The statement of Theorem II .11 can be reformulated 
as follows: H n+ \{x, y) = if and only if W £ Cj~ for some k < n. 

Acknowledgments. I would like to thank Alberto Griinbaum for his 
collaboration in jlUj which inspired my interest in the subject. I also want 
to thank George Wilson for suggesting this simple proof of the implication 
(i)=Kii). 



References 

[I] M. Adler and J. Moser, On a class of polynomials connected with the Korteweg-de 
Vries equation, Comm. Math. Phys. 61 (1978), 1-30. 

[2] H. Airault, H. P. McKean and J. Moser, Rational and elliptic solutions of the Korteweg- 
de Vries equation and a related many-body problem, Comm. Pure Appl. Math. 30 
(1977), 95-148. 

[3] I. Avramidi and R. Schimming, A new explicit expression for the Korteweg-de Vries 
hierarchy, Math. Nachr. 219 (2000), 45-64. 

[4] Y. Berest, Huygens principle and the bispectral problem, In: The bispectral problem 
(Montreal, PQ, 1997), 11-30, CRM Proc. Lecture Notes, 14, Amer. Math. Soc, Prov- 
idence, RI, 1998. 

[5] Y. Berest and A. P. Veselov, The Huygens principle and integrability, (Russian) Us- 
pekhi Mat. Nauk 49 (1994), 7-78, translation in Russian Math. Surveys 49 (1994), 
5-77. 

[6] J. Burchnall and T. Chaundy, Commutative rings of ordinary differential operators. 

II. The identity P" = Q m , Proc. Royal Soc. London (A) 134 (1931), 471-485. 
[7] E. Date, M. Jimbo, M. Kashiwara and T. Miwa, Transformation groups for soliton 

equations, in: M. Jimbo and T. Miwa (eds), Proc. RIMS Symp. Nonlinear Integrable 

Systems - Classical and Quantum Theory (Kyoto 1981), World Scientific, Singapore 

1983, pp. 39-119. 

[8] J. J. Duistermaat and F. A. Griinbaum, Differential equations in the spectral parame- 
ter, Comm. Math. Phys. 103 (1986), 177-240. 

[9] F. A. Griinbaum, Some bispectral musings, In: The bispectral problem (Montreal, PQ, 
1997), 11-30, CRM Proc. Lecture Notes, 14, Amer. Math. Soc, Providence, RI, 1998. 

[10] F. A. Griinbaum and P. Iliev, Heat kernel expansions on the integers, Math. Phys. 
Anal. Geom. 5 (2002), 183-200, |math.CA/ 0206089 

[II] J. Hadamard, Lectures on Cauchy's Problem, New Haven, Yale Univ. Press, 1923. 



FINITE HEAT KERNEL EXPANSIONS ON THE REAL LINE 



7 



[12] L. Haine, The spectral matrices of Toda solitons and the fundamental solution of some 
discrete heat equations, to appear in Annales de Hnstitut Fourier. 

[13] R. Hirota, The direct method in soliton theory, Translated from the 1992 Japanese 
original and edited by Atsushi Nagai, Jon Nimmo and Claire Gilson (with a foreword 
by Jarmo Hietarinta and Nimmo) . Cambridge Tracts in Mathematics, 155. Cambridge 
University Press, Cambridge, 2004. 

[14] P. Iliev, On the heat kernel and the Korteweg-de Vnes hierarchy, to appear in Annales 
de l'lnstitut Fourier, available at math-ph/0504045 

[15] M. Kac, Can one hear the shape of a drum?, Amer. Math. Monthly 73 (1966), 1-23. 

[16] N. Lebedev, Special functions & their applications, Dover Publications, Inc., 1972. 

[17] H. P. McKean and P. van Moerbeke, The spectrum of Hill's equation, Invent. Math. 
30 (1975), 217-274. 

[18] M. Sato and Y. Sato, Soliton equations as dynamical systems on infinite dimensional 
Grassmann manifolds, Lect. Notes Num. Appl. Anal. 5 (1982), 259-271. 

[19] R. Schimming, An explicit expression for the Korteweg-de Vries hierarchy, Acta Appl. 
Math. 39 (1995), 489-505. 

[20] G. Segal and G. Wilson, Loop groups and equations of KdV type, Publ. Math. IHES 
61 (1985), 5-65. 

School of Mathematics, Georgia Institute of Technology, Atlanta, GA 
30332-0160, USA 

E-mail address: iliev@math.gatech.edu 



